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■ Abstract 

A fully differential calculation of the next-to-leading order QCD corrections to the production of 
Z-boson pairs in association with a hard jet at the Tevatron and LHC is presented. This process 
is an important background for Higgs particle and new physics searches at hadron colliders. We 

ON 

find sizable corrections for cross sections and differential distributions, particularly at the LHC. 



Residual scale uncertainties are typically at the 10% level and can be further reduced by applying 



a veto against the emission of a second hard jet. Our results confirm that NLO corrections do not 
simply rescale LO predictions. 

PACS numbers: 12.38.Bx, 13.85.-t, 14.70.Hp 
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I. INTRODUCTION 



Weak boson pair production at hadron colliders plays an essential part in the search for 
Higgs particles and for physics beyond the Standard Model (SM), since weak bosons can 
decay into jets, charged leptons or neutrinos and hence produce the same signatures as Higgs 
bosons, new coloured particles, new electroweak gauge bosons or dark matter candidates. 
In addition to being an important background to direct new physics searches at the Large 
Hadron Collider (LHC) |l( > weak boson pair production also allows to search for new physics 
via experimental evidence for SM deviations in the form of anomalous interactions between 
electroweak gauge bosons [2J]. 

Z-boson pair production has been observed at the Tevatron [3| and studied extensively 



by the LHC general-purpose detector collaborations [4J]. Since LO predictions for hadron 
collider processes are affected by large QCD scale uncertainties with respect to normalisation 
and kinematical dependence, the inclusion of NLO QCD corrections is important when 
comparing predictions for cross sections and differential distributions with data. Theoretical 
predictions for ZZ production at leading order (LO) \^ have thus been improved by including 
next-to-leading order (NLO) QCD corrections without [6] and with decays 7]. More recently, 
Z-boson pair production at NLO has also been investigated in selected SM extensions 8]. 
As a window to new physics ZZ production is particularly interesting because of the absence 
of ZZ^i and ZZZ couplings [9] in the SM. Probing such anomalous neutral gauge boson 



couplings at hadron colliders has also been studied in the literature [10]. 

Going beyond final states with two particles, NLO QCD corrections have been computed 



for the production of three wea 
association with up to three jets 



c/yector bosons [ll], the production of a weak boson in 



12j l and weak boson pair production in vector boson fusion 



[131 ] . Of particular interest is also the production of weak boson pairs with one additional 
jet at NLO. This process is interesting in its own right, due to the enhanced jet activity, 
particularly at the LHC. It also provides the real- virtual contribution to the next-to- next- 
to-leading order (NNLO) corrections to weak boson pair production. The production of 
H^-boson pairs with an additional jet has thus been calculated at NLO without [LjJ and 



with 



15 



16| decays. An additional contribution to the NNLO corrections that has been 



calculated for WW and ZZ production (at LO) including decays is the loop-induced gluon- 
fusion process 17j. Other building blocks for the NNLO calculation of weak boson pair 



2 



fWWW z Q 



•(WWW z 



fWWW z 



'^RRTOTT g 



"0151T0TP g q 



•WWW Z 



FIG. 1: Representative LO graphs for the partonic process qq — ► Z-Zg. 

production have been presented in Ref. 

In this paper, we present a calculation of the 0(a s ) corrections to Z-boson pair production 
with an additional jet at hadron colliders in the SM. 1 Details of the NLO calculation are 
described in Sec. [Til We then present numerical results in Sec. IIIII and end with conclusions. 



II. DETAILS OF THE NLO CALCULATION 



At LO, all channels for ZZj production at hadron colliders are related to the amplitude 
— > ZZqqg by crossing symmetry. Therefore, the following subprocesses contribute: 

qq -> ZZg , qg -> ZZq , qg -> ZZq , 

where q can be either an up- or down-type quark. 2 We calculate in the 5-flavour scheme, 
i.e. q = u,c, d, s, b, and neglect all quark masses. Representative LO diagrams for the first 
subprocess are shown in Fig. [TJ Comparing the ZZ with the corresponding WW production 
amplitude, key differences are that the H^-boson coupling to fermions is purely left-handed, 
which leads to a reduced number of helicity amplitudes in that case, the produced weak 
bosons are distinct for WW production, but identical for ZZ production (leading to a 
significant increase in the number of Feynman diagrams), and graphs with triple-gauge 
boson vertices contribute to WW, but not to ZZ production. 



1 Partial results of our calculation have already been presented in Ref. 19 1. 

2 The down-type quark initiated amplitudes are obtained from the up-type quark initiated amplitudes by 
adjusting the chiral couplings. 
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FIG. 2: Representative one-loop graphs for the partonic process qq — ► Z-Zg. 



Virtual corrections 



At (D(a s ), the most complicated loop topologies are pentagon graphs derived from the 
tree-level graphs via virtual gluon exchange (and crossing) and box graphs derived by closing 
the quark line in the tree- level graphs and attaching a gqq current. Since we calculate with 
Nf = 5 and neglect quark mass effects, graphs with Higgs boson exchange do not contribute. 
Triangle graphs, where three gauge bosons couple to a quark loop, also vanish in the massless 
quark limit. Representative one- loop graphs for the partonic process qq — > ZZg are shown 
in Fig. [2j Starting from the Feynman graph representation, two independent sets of ampli- 
tude expressions have been generated: one manually, the other using QGRAF [20(]. Both 
representations employ the spinor helicity formalism of Ref . 21] . Polarisation vectors have 
been represented via spinor traces, i.e. kinematic invariants up to global phases. By obtain- 
ing an analytical representation for the full amplitude, we aim at promoting simplification 
via analytical cancellations. Especially we employ that, apart from the rank one case, all 
pentagon tensor integrals are reducible, i.e. can directly be written as simple combinations 
of box tensor integrals. For the remaining tensor integrals we employ the GOLEM-approach 



22j . In this approach, the use of 6-dimensional IR finite box functions allows to isolate IR 



divergences in 3-point functions. We use FORM [23] and Maple to obtain tractable analyti- 
cal expressions for the coefficients to the employed set of basis functions for each independent 
helicity amplitude, and to further simplify them. The basis functions are evaluated using 
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the GOLEM95 implementation 2jj. We note that for the reduction of box topologies one 



25|. If one fully 



obtains the same result as with the Passarino-Veltman tensor reduction 
reduces all tensor integrals to a scalar integral representation, the difference between the 
two approaches results from the treatment of the pentagon integrals and the use of finite 
6-dimensional box functions. 



To treat 75 we employ the 't Hooft-Veltman scheme 26], where the 7^ are split into a 
4-dimensional part that anti-commutes with 75 and a commuting remainder. 3 As is well 
known, to take into account differences between the QCD corrections to axial vector and 
vector currents, a finite renormalisation has to be performed. To enforce the correct chiral 
structure of the amplitudes, a finite counterterm for the axial part is included in the used 
gauge boson vertex (see e.g. Ref. [281]): 



Vv q q -9vl' 1 + Z 5 g a 7^75 with Z 5 = 1 - C F 



a. 



7T 



We have verified that the relative contribution of graphs with quark loops to the inte- 
grated results is well below 1%. We therefore neglect this contribution in the following. 
We have compared our two independent implementations of the virtual amplitudes - both 
generated using the GOLEM reduction - and have found agreement of 9-16 significant digits 
for all contributions at two test phase space points. The discrete Bose, P and C symmetries 
induce relations between different helicity amplitudes that have been verified numerically 
as additional check. We have furthermore tested gauge invariance by confirming the Ward 
identity for external gluons. Furthermore, we used the same tools as in the calculation pre- 
sented here to calculate numerical results for the successful comparison of virtual corrections 
to WWj production in Ref. 29]. That comparison therefore provides an additional check 
of our calculation. We have also verified that our LO amplitude implementation is correct. 
To calculate numerical results for the virtual contributions we employed the OmniComp- 
Dvegas package [30j , which facilitates parallelised adaptive Monte Carlo integration and was 



developed in the context of Ref. 



3l|. 



3 Note that the 't Hooft-Veltman scheme treats the observed particles in 4 dimensions, but the soft/collinear 
gluons in d dimensions. This guarantees that for the IR subtractions the same Catani-Seymour dipole 
terms as for conventional dimensional regularisation can be used [27| . 
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B. Real corrections 



The 0(a s ) real corrections channels for ZZj production at hadron colliders are related 
to the amplitudes — > ZZqqgg and — > ZZqqq'q' by crossing symmetry. While all virtual 
corrections channels are already active at LO, new real corrections channels open up at 
NLO, namely the gg, qq', qq' (q' 7^ q) and qq' channels. Note that these new channels are 
effectively of LO type. 

To facilitate the cancellation of soft and collinear singularities we employ the Catani- 



Seymour dipole subtraction method [32]. We use the SHERPA implementation 33J to 
calculate numerical results for the finite real corrections contribution. All amplitude and 
dipole contributions have been verified through comparison with results calculated with two 



independent implementations: 4 MadDipole/MadGraph 351] and HELAC 36 1. In both com- 
parisons 9-significant-digits agreement or better was achieved for all contributions for two 
test phase space points. We have also successfully compared with an in-house implementa- 
tion of a set of independent dipoles. 

III. NUMERICAL RESULTS 

In this section, we present first NLO predictions for ZZj cross sections and differential 
distributions at the Tevatron (pp, y/s = 1.96 TeV) and LHC (yy, J~s = 14 TeV) and compare 



5 



them with the NLO results for WWj production given in Ref. [la].'- 

As mentioned above, our calculation employs the 5-flavour scheme and the massless quark 
approximation (including the b quark). We use the SM parameters 6 

M z = 91.188 GeV, a(M z ) = 0.00755391226 , sin Q 2 W = 0.222247, 



and employ CTEQ6 parton distribution function sets [37j]. LO (NLO) cross sections are 
calculated with CTEQ6L1 (CTEQ6M) and LO (NLO) a s running. For a s (M z ), the default 
LHAPDF H values are used: a fl (91.188 GeV) = 0.129783 at LO and a s (91.70 GeV) = 
0.1179 at NLO. In our parton-level calculation, partons are clustered into jets using the 
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34] 



Other implementations that automate the dipole subtraction method of Ref. [321 ] have been reported in 
Ref. 



5 Differences in the input parameters are minute (less than 0.1%). 

6 We provide the exact input values of our calculation in order to facilitate reproducibility. 
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FIG. 3: Comparison of the scale dependence (fiR = f^F = A*) of the ZZ+jet cross section at the 
Tevatron and LHC with prjet > 50 GeV for the hardest jet in LO (dotted) and NLO (solid). For 
the LHC the exclusive NLO cross section when a £>T,jet > 50 GeV veto for additional jets is applied 
is also shown (dot-dashed). Input parameters are defined in the main text. 



inclusive k t algorithm 39] with R = 0.7. To study the scale dependence of cross section pre- 
dictions we use the discrete grid /i = 2 % l 2 Mz with i E {—7, —6, . . . , 7}. The renormalisation 
and factorisation scales are identified here (/i = fiR = fip)- We apply a > 50 GeV cut on 
the hardest jet unless noted otherwise. 

In Fig. [31 LO and NLO predictions for ZZj production cross sections at the Tevatron 
and LHC are displayed as function of the QCD scale /i, which is varied by a factor 10 around 
the Z mass. At LO, we find a much larger scale variation at the Tevatron than at the LHC. 
When NLO corrections are included, the Tevatron cross section reaches its maximum at 
approximately Mz/2 and its variation is very effectively reduced. The shape of the cross 
section variation at the LHC on the other hand is qualitatively unchanged when going from 
LO to NLO. We attribute this to new channels that become active at NLO (see Sec. Ill Bh . 
which have a modest impact at the Tevatron, but a sizable impact at the LHC, due to 
parton densities being probed in different x regions. One might be tempted to conclude 
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Aa/a(pp - 


-> ZZ + jet), y/s -- 


= 1.96 TeV 




H/Mz G [|,2] 


M/Mz6[l,4] 


/x/M z e[|,8] 


LO 


23% 


44% 


62% 


NLO 


6% 


11% 


19% 



TABLE I: Scale uncertainty for LO and NLO cross sections for ZZ+jet production at the Tevatron 
as function of the scale variation. The relative scale uncertainty is defined through the envelope: 
Aa/o := [cr ma _ x (p £ I) - a min (fi G I)]/[o- max (/i <E /) + fx min (/i £ I)}. Other details as in Fig. El 





Aa/a(pp 


-> ZZ + jet), y/s 


= 14 TeV 




M/MzG[i,2] 


M/MzG[i,4] 


n/M z e[l,8] 


LO 


12% 


23% 


34% 


NLO 


7% 


15% 


23% 


NLO with 2 nd jet veto 


0.5% 


3% 


6% 



TABLE II: Scale uncertainty for LO and NLO cross sections for ZZ + jet production at the LHC 
as function of the scale variation. Other details as in Fig. [3] and Tabled 



that a constant If -factor is a good approximation for the inclusive NLO cross section at 
the LHC. However, the K factor (also shown) varies between 1.3 and 1.6 in the displayed 
scale range. Following Ref. [lij], we also calculate an exclusive NLO cross section for the 
LHC by vetoing 2-jet events with a second hardest jet with > 50 GeV (NLO with 2 nd 
jet veto). This exclusive NLO LHC cross section decreases for scales below Mz and has a 
strongly reduced scale uncertainty. Since the qualitative features of our results are similar 
to those found for WWj production, we quantify the ZZj LO and NLO scale uncertainties 
for the Tevatron and LHC in Tables d and [Til respectively. When comparing our results 
for fi/Mz G [|,2] with the corresponding WWj results extracted from Tables 4 and 1 in 



Ref. [16| , we find deviations of less than 2 percentage points. 

In Tables IIHI and IIVI we show for the Tevatron and LHC, respectively, how the ZZj 
LO and NLO cross sections change when the pt cut on the hardest jet is varied. These 
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a(pp — > ZZ + jet) [pb], y/s = 1.96 TeV 


Prjet cut [GeV] 


20 


50 


100 


200 


LO 
NLO 


0.27202(3) 
0.3307(6) 


0.07456(1)1^ 
0.0836(1)^5% 


0.016037(2) 
0.01583(4) 


0.0012651(1) 
0.000976(4) 



TABLE III: ZZ+jet production cross section at the Tevatron with different pr cuts for the hardest 
jet. The scale /i = Mz is employed with [ir = fip = \i. The integration error is given in brackets. 
The minimum and maximum relative deviation from a(n = Mz) obtained with the scale variation 
n/Mz G [|)2] is shown as sub- and superscript for a prjet cut of 50 GeV. Input parameters are 
defined in the main text. 





a(pp -> ZZ + jet) [pb], = 14 TeV 


PTjet cut [GeV] 


20 


50 


100 


200 


LO 
NLO 

NLO with 2 nd jet veto 


6.505(1) 
8.01(3) 


2.6978(4)1^ 
3.653(9)lg 

2.637(9)15% % 


1.0066(1) 
1.511(4) 
0.755(4) 


0.22974(3) 
0.415(2) 
0.1005(9) 



TABLE IV: ZZ+jet production cross section at the LHC with different pt cuts for the hardest 
jet. For cut values above 20 GeV, we also give the NLO cross section when a PTjet > 50 GeV veto 
for additional jets is added to the selection. Other details as in Table IIIII 

results demonstrate that all K factors are pt cut dependent. In general, the K factor for 
ZZj production will have a non-negligible differential dependence. As example, we display 
in Fig. H]the differential LO and NLO distributions with respect to the invariant ZZ mass 
and the resulting K factor at the Tevatron and LHC. The i^-factor bands shown in this 
figure correspond to a variation of the scale /x by a factor of 2 in the NLO differential cross 
section only, i.e. we display [da^j J o/dMzz](f^)/[dcri,o/dMzz}(Mz) with fi/Mz G [§,2]. One 
can distinguish the modest variation of the inclusive NLO K factor at the LHC from the 
strong decrease of the other K factors with increasing M Z z- 
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FIG. 4: ZZ invariant mass distribution for ZZ+jet production at the Tevatron and LHC with 
= fJ'F = Mz- The differential K factor is also shown. The iT-factor bands are defined in the 
main text. Other details as in Fig. [3j 



IV. CONCLUSIONS 



In this paper we have presented first results for ZZj production at NLO QCD accuracy, 
obtained with fully differential parton-level Monte Carlo programs that allow to take into 
account realistic experimental selection cuts. For a default scale choice of /x = Mz and a 
Pt cut of 50 GeV for the hardest jet we find a K factor of 1.1 and 1.35 at the Tevatron 
and LHC, respectively. At the Tevatron, the NLO corrections stabilise the LO prediction 
for cross sections considerably. The shape of the cross section variation at the LHC on the 
other hand is qualitatively unchanged when going from LO to NLO. Nevertheless, at the 
LHC stabilisation at NLO can still be achieved by applying suitable selection cuts like for 
example a veto against the emission of a second hard jet. Our results indicate that residual 
scale uncertainties are typically at the 10% level and can be further reduced to about 5% or 
less by applying suitable selections. 
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